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Abstract 

We investigate the circular motion of charged test particles in the gravitational field of a charged 
mass described by the Reissner-Nordstrom (RN) spacetime. We study in detail all the spatial 
regions where circular motion is allowed around either black holes or naked singularities. The 
effects of repulsive gravity are discussed by finding all the circles at which a particle can have 
vanishing angular momentum. We show that the geometric structure of stable accretion disks, 
made of only test particles moving along circular orbits around the central body, allows us to 
clearly distinguish between black holes and naked singularities. 
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I. INTRODUCTION 



Let us consider the background of a static gravitational source of mass M and charge Q, 
described by the Reissner-Nordstrom (RN) hne element in standard spherical coordinates 

ds^ = -^de + ^cir^ + {de^ + sin^ Odcf)'') , (1) 

where A = (r - r+)(r - r_) and r± = M ± M"^ - Q2 

are the radii of the outer and inner 
horizon, respectively. Furthermore, the associated electromagnetic potential and field are 

A = ^dt, F = dA = -^dtAdr , (2) 

respectively. 

The motion of a test particle of charge q and mass fi moving in a RN background ([T]) is 
described by the following Lagrangian density: 

^ = ^ga^i'^i'' + eAo^x", (3) 

where are the components of the electromagnetic 4-potential, the dot represents differ- 
entiation with respect to the proper time, and the parameter e = q/fi is the specific charge 
of the test particle. The equations of motion of the test particle can be derived from Eq. ([3]) 
by using the Euler-Lagrange equation. Then, 

i;"VaX^ = eF'^^x\ (4) 

where F^p = A^^^p - Ap^^. 

Since the Lagrangian density ([3]) does not depend explicitly on the variables t and 0, the 
following two conserved quantities exist 



dC /A. eQ\ E 

dt ^ / 



= ^ = r^ sin^ 9(1) = -, (6) 
d(p 

where L and E are respectively the angular momentum and energy of the particle as mea- 
sured by an observer at rest at infinity. Moreover, to study the motion of charged test 
particles in the RN spacetime it is convenient to use the fact if the initial position and the 
tangent vector of the trajectory of the particle lie on a plane that contains the center of the 
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body, then the entire trajectory must he on this plane. Without loss of generality we may 
therefore restrict ourselves to the study of equatorial trajectories with 9 = 7r/2. 

On the equatorial plane 6 = 7r/2, the motion equations can be reduced to the form 
j.2 _^ Y2 _ 1 ^2 describes the motion inside an effective potential V. Then, we 

define the potential 

/i r V \ n'^r'^ J \ r J 

as the value oi E / ^ that makes r into a "turning point" {V = E/fi); in other words, the 
value of at which the (radial) kinetic energy of the particle vanishes jl- J]. The effective 
potential with positive (negative) sign corresponds to the solution with 

lim = ( lim E^ = — /i ) , 

where 

E+{L,e,r)>E-{L,e,r), (8) 

and the following relation holds: 

E^iL,e,r) = -E'iL,-e,r). (9) 
The behavior of the effective potential strongly depends on the sign of eQ; in particular in 



the case of < 0, negative energy states for the solution E~^ can exist (see also |5l-lll|). 



The problem of finding exact solutions of the motion equations of test particles moving 
in a RN spacetime has been widely studied jn literature in many contexts and ways. For a 



recent discussion we mention the works |5Nll|. In particular, in a recent paper |5| the full 



vays 



set of analytical solutions of the motion equations for electrically and magnetically charged 
test particles is discussed in terms of the Weierstrass (7, a and () functions. The general 
structure of the geodesies was discussed and a classification of their types was proposed. 
Remarkably, analytical solutions are found in the case of a central RN source not only with 
constant electric charge, but also with constant magnetic charge. It is interesting to notice 
that if either the test particle or the central body possesses both types of charge, it turns 
out that the motion is no longer confined to a plane. In the present work, we consider 
only equatorial circular orbits around a central RN source with constant electric charge. 
Instead of solving directly the equations of motion, we explore the properties of the effective 
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potential function associated to the motion. Thus, we discuss and propose a classification of 
the equatorial orbits in terms of the two constants of motion: the energy E/ ^ and the orbital 
angular momentum L/{^M). In fact, we focus our attention on some peculiar features of the 
circular motion and the physics around black holes and naked singularities. In particular, 
we are interested in exploring the possibility of distinguishing between black holes and naked 
singularities by studying the motion of circular test particles. In this sense, the present work 



complements and is different from previous studies 
In a previous work 



M- 



1^ . Q , we analyzed the dynamics of the RN spacetime by studying 
the motion of neutral test particles for which the effective potential turns out to coincide with 
as given in Eq.([7]) with e = 0. We will see that in the case of charged test particles the 
term eQ/r drastically changes the behavior of the effective potential, and leads to several 
possibilities which must analyzed in the case of black holes and naked singularities. In 
particular, we will show that for particles moving along circular orbits there exist stability 
regions whose geometric structure clearly distinguishes naked singularities from black holes 
(see also Q, Q and [igI, The plan of this paper is the following: In Sec. [Ill we 

investigate the behavior of the effective potential and the conditions for the motion of positive 
and negative charged test particles moving on circular orbits around the central charged 
mass. This section also contains a brief analysis of the Coulomb approximation of the 
effective potential. In Sec. IIIIl we will consider the black hole case while in Sec. [IV]we shall 
focus on the motion around naked singularities. The conclusions are in Sec. |Vl 

II. CIRCULAR MOTION 

The circular motion of charged test particles is governed by the behavior of the effective 
potential (171). In this work, we will mainly consider the special case of a positive solution 
V+ for the potential in order to be able to compare our results with those obtained in the 



case of neutral test particles analyzed in 
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131]. Thus, the radius of circular orbits and 



the corresponding values of the energy E and the angular momentum L are given by the 
extrema of the function Therefore, the conditions for the occurrence of circular orbits 
are: 

^^0. V,^^. (10) 
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When possible, to simplify the notation we will drop the subindex (+) so that, for example, 
V = E/fi will denote the positive effective potential solution. Solving Eq. (|T0|) with respect 
to L, we find the specific angular momentum 

^ = ^ [2(Mr - g2)S + e'Q'A ± QA^e^ (4S + e^g^)] , (n) 

where S = — 3Mr + 2Q^, of the test particle on a circular orbit of radius r. The corre- 
sponding energy can be obtained by introducing the expression for the angular momentum 
into Eq. ([7]). Then, we obtain 



E^_eQ A^2i: + e^Q^±Q^e\AT. + e^Q^) 
~ r V2r|S| 
The sign in front of the square root should be chosen in accordance with the physical 
situation. This point will be clarified below by using the formalism of orthonormal frames. 

An interesting particular orbit is the one in which the particle is located at rest as seen 
by an observer at infinity, i.e., L = 0. These "orbits" are therefore characterized by the 
following conditions 

1 = 0. f = 0. (13) 



18| . Solving Eq. ( lT3l) for Q ^ and e 7^ 0, we find the following radius 



Table [I shows the explicit values of all possible radii for different values of the ratio Q/M. 
A particle located at r = with angular momentum L = will have the energy (see also 
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fl g 1 V e2 - 1 , / e2(Af2-Q2)(,2_-T) - I " ^ ^ 



The minimum radius for a stable circular orbit occurs at the inflection points of the 
effective potential function; thus, we must solve the equation 



0, (16) 



for the orbit radius r, using the expression f lTT]) for the angular momentum L. From Eq. f lTOj) 
and Eq. ( IT6|) we find that the radius of the last stable circular orbit and the angular mo- 
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TABLE I: Radii of the "orbits" characterized by the conditions L = and dV/dr ~ 0. 



mentum of this orbit are related by the following equations 



2\ ^3 



and 



where in order to simplify the notation we introduced the normalized quantities L — )■ 
L/{M/fi), r — r/M, and Q — t- Q/M. Equation ( ITTl) depends on the test particle spe- 
cific charge e via the function L as given in Eq. ( ITTj) . It is possible to solve Eq. ( IT7|) for the 
last stable circular orbit radius as a function of the free parameter L. We find the expression 

(Lf )2 r2 

V Isco/ 



[2g^(5M - 3r)r - SQ^ - r2[6M2 + (r - 6M)7 



/x2 2 [2g4 + 3gV(r - 2M) - (2r - 3M)r2] 

±\/9Q2 + (r - 6M)r (Q^ + (r - 2M)r)^^^ 



for the angular momentum of last stable circular orbit. Eq. ( IT8|) can be substituted in 
Eq. f|T7|) to find the radius of the last stable circular orbit. 



A. Coulomb potential approximation 

Consider the case of a charged particle moving in the Coulomb potential 
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This means that we are considering the motion described by the following effective potential 

E+ eQ I , , 

^+ = — = — + 19 

where eQ < 0. The Coulomb approximation is interesting for our further analysis because 
it corresponds to the limiting case for large values of the radial coordinate r [cf. Eq.([7])]. 
Circular orbits are therefore situated at r = with 




and — > e^Q^ (20) 

and in the case e = with Q > 0, circular orbits exist in all r > for L = 0. We conclude 
that in this approximation circular orbits always exist with orbital radius Tc and angular 
momentum satisfying the condition > |eQ|- For the last stable circular orbit situated 

at r = rigj.o we find 

= with ?L^^J^ = and = \eQ\ . (21) 

/i yU 

This means that, in the approximation of the Coulomb potential, all the circular orbits are 
stable, including the limiting case of a particle at rest on the origin of coordinates. 

Furthermore, Eqs. ( l20] - [2Tl) show that, in contrast with the general RN case, for a charged 
particle moving in a Coulomb potential only positive or null energy solutions can exist. See 
Fig. [1] where the potential ( lT9l) is plotted as a function of the orbital radius for different 
values of the angular momentum. 



III. BLACK HOLES 



In the case of a black hole (M^ > Q^) the two roots V± of the effective potential are 
plotted as a function of the ratio r/M in Fig. |2] for a fixed value of the charge-to-mass 
ratio of the test particle and different values of the angular momentum L/(M/i) (see also 



25 



30|). Notice the presence of negative energy states for the positive solution V+ = j \i 
of the effective potential function. Negative energy states for V+ are possible only in the 
case eQ < 0. In particular, the largest region in which the V+ solution has negative energy 
states is 



M + ^/M^ - g2 < r < M + ^/W^^Q^T^^ (22) 
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FIG. 1: Effective potential for a cliarged test particle with e = —2 moving in a Coulomb potential with 
Q/M = 2 for different values of the momentum L* = L/{iiM). The points indicate the minima of the 
potential. In particular, for L* = \eQ\/M the potential vanishes on the origin r = (see text). 

and corresponds to the limiting case of vanishing angular momentum (L = 0). For L 7^ 
this region becomes smaller and decreases as L increases. For a given value of the orbit 
radius, say tq, such that ro < M + Y/M2~ir"Q2"(YTr"^^ i\^q angular momentum of the test 
particle must be chosen within the interval 



for a region with negative energy states to exist. This behavior is illustrated in Fig. [2l 

Fig. [3] shows the positive solution V+ of the effective potential for different values of the 
momentum and for positive and negative charged particles. In particular, we note that, at 
fixed Q/M for a particle with |e| < 1, in the case eQ > the stable orbit radius is larger 
than in the case of attractive electromagnetic interaction, i. e., eQ < 0. In Fig. HI the 
potential V+ of an extreme black hole is plotted for different, positive and negative values 
of the test particle with charge-to-mass ratio e. In this case, it is clear that the magnitude 
of the energy increases as the magnitude of the specific charge of the particle e increases. 

As mentioned in Sec. Ull in the case of the positive solution for the effective potential the 
conditions for the existence of circular orbits 




(23) 



0, V 



E dV 



(24) 



r = 



11 ' dr 
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FIG. 2: The effective potential as a function of r/M for a charged particle of charge-to-mass ratio e = q/^ 
moving in a Reissner-Nordstrom black hole of charge Q and mass M . The graphic shows the positive E'^ j [i 
(black curves) and negative roots / ji (gray curves) of the effective potential for Q/il/ = 0.5, e = —2, and 



different values of the momentum L* = L/(MiJ.). The outer horizon is located at r+ = M + \J A/^ — ~ 
1.87Af . Note the presence of negative energy states for the positive roots. 

lead to Eqs. ffTT]) and f|T2l) in which the selection of the (±) sign inside the square root should 
be done properly. To clarify this point we consider explicitly the equation of motion for a 
charged particle in the gravitational field of a RN black hole. 

a(f/)° = eF'^f^U^ (25) 

where a(t/) = Vf/t/ is the particle's 4-acceleration. Introducing the orthonormal frame 



Ai/2 



dr, e-Q = -dg, 67 = —— 
r J. V J. gj]^ y 



with dual 



A^/^ - r ' 

* — dt, u"" = ^^dr, = rdO, u'^ = r sin 



r Ai/2 
the tangent to a (timelike) spatially circular orbit u° can be expressed as 



(26) 



(27) 



u = T{dt + (d^) = 7 (e£ + z^e^) 



where F and 7 are normalization factors 



T' = {-gu-e9Hr' and f = (1 - i.^)-!^ 
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FIG. 3: The effective potential for a charged particle of charge-to-mass ratio, e — q/fi, moving in a 
Reissner- Nordstrom spacetime of charge Q and mass M with charge-to-mass ratio Q/M = 0.5 is plotted as 
a function of the radial coordinate r/M for different values of the angular momentum L* = L/{AIfi). The 
outer horizon is located at r+ w 1.87A/. In the graphic on the left with e = 0.1, the effective potential for 
L* « 3.5 has a minimum Vmin ~ 0.954 at r„ii„ « 8.84A'/. In the graphic on the right with e — —0.1, the 
minimum Vmin ~ 0.94 is located at r„im ~ 7.13M for L* « 3.5. 
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FIG. 4: The effective potential V+ is plotted as a function of r/M for a charged test particle with specific 
charge e = q/ ji moving in the field of a Reissner- Nordstrom extreme black hole {Q ~ M). Here L/{Mfj,) = 4, 
and the effective potential is plotted for different values of e. The outer horizon is located at r+ = M + 
\J — (52 — Note the presence of negative energy states for particles with negative e. 

which guarantees that Uau"- = — 1. Here C is the angular velocity with respect to infinity 
and V is the "local proper linear velocity" as measured by an observer associated with the 
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orthonormal frame. The angular velocity ( is related to the local proper linear velocity by 



> / 9tt 



Since only the radial component of the 4-velocity is non-vanishing, Eq. f l2^ can be written 
explicitly as 



Ug eQ 



where 



c 



Mr -Q2 



9 ^ ^2 ^ '^9 



Mr-Q^ 



A 



(28) 



(29) 



This equation gives the values of the particle linear velocity u = '^vf which are compatible 
with a given value of eQ on a circular orbit of radius r, i. e., 



± 



\ 2t%1 r\gVg y 72 ' ^r%l 



(30) 



where 



and 



I9 



r2 - 3Mr + 2g2 



1/2 



it = (1 - ^, 



±2\-l/2 



(31) 



In the limiting case of a neutral particle (e = 0), Eq. (l28|l implies that the linear velocity of 
the particle is Vg. 

We introduce the limiting value of the parameter e corresponding to a particle at rest, 
= 0, in Eq. ([2H]), i. e.. 



Mr - Q' 



By introducing this quantity into Eq. fl28|) . one gets the following equivalent relation 

- = 7 1 - — 
whose solution (1301] can be conveniently rewritten as 



(32) 



(33) 



A ± VA2 - 1 + (e/ 



1/2 



where 



2 Veo 



(34) 



(35) 
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Moreover, from Eq. (l33l) it follows that e < implies that > (because eo is always 
positive for r > r+), so that the allowed solutions for u can exist only for r > r+, where 



r 

7 



= ^(3M+ v/9M2 _ 8Q2 j ^ (36) 



the equality corresponding to Ug = 1. In this case, the solutions of Eq. fl28|) are given by 
u = ±z/+. 

For e > 0, instead, solutions can exist also for < r < . The situation strongly 
depends on the considered range of values of e and is summarized below. 
Equation (1311) gives the following conditions for the existence of velocities 

- 1 + {e/eof > , (37) 
A ± ^A2 - 1 + (e/eo)2 > . (38) 

The second condition, Eq. (1551) . is satisfied by 

3M 1 



r>ri = — + -v/9M2-8Q2_e2Q2_ (39) 

Moreover for Q = M and e = 1 it is A + a/A^ - 1 + (e/eo)^ > when M < r < (3/2)M. 
However it is also possible to show that condition Eq. ( TTOl) is satisfied for < Q < M and 
e > only in the range r > r^. 

Requiring that the argument of the square root be nonnegative implies 

.<,.^/^MEW. (40) 

The condition ( !38l) will be discussed later. 

From the equation of motion (l33i) it follows that the velocity vanishes for e/eo = 1, i. e., 
for [cf. Eq.dTID] 



r = r. 



e2g2 _ ^2 [ 

which exists only for e > M/Q. We thus have that 

e 



M(e2-1) + v/e2(e2-l)(M2-Q2) ^ (41) 



whereas 



> 1 for r > r^, (42) 
Co 



— < 1 for r+ < r < . (43) 
Co 
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On the other hand, the condition u = in Eq. (IMI) imphes that 



-I e/eo = l 



1. e. 



0, 
0, 



thus z/^ is identically zero whereas 

1 



V2Q 



- 2A{rs) = only for 
5M2 - 4Q2 + ^25M2 - 24g2. 



(44) 
(45) 



(46) 



Finally, the lightlike condition z/ = 1 is reached only at r = r+, where = 1 = u. 

The behavior of charged test particles depends very strongly on their location with respect 
to the special radii r_|_, ri, r+, and r^. In Sec. IIII Al the behavior of these radii will be analyzed 
in connection with the problem of stability of circular orbits. 

On the other hand, the particle's 4-momentum is given by P = mil — qA. Then, the 
conserved quantities associated with the temporal and azimuthal Killing vectors ^ = dt and 
Tj = ds are respectively 



r 

r 

P ■ ri = — 7Z/ 



1- 



r 
L 



E 



(47) 
(48) 



where E/ ^ and L//i are the particle's energy and angular momentum per unit mass, respec- 
tively (see also Eqs.(III]) and f|T2l) ). 
Let us summarize the results. 



1. Case e < 

The solutions are the geodesic velocities v = in the range r > r+ as illustrated 
in Fig. |5l Orbits with radius r = r+ are lightlike. We can also compare the velocity of 
charged test particles with the geodesic velocity Ug for neutral particles. For r > we see 
that > Vg always. This means that, at fixed orbital radius, charged test particles acquire 
a larger orbital velocity compared to that of neutral test particles in the same orbit. As it 
is possible to see from Eq. ( |30l) and also in Fig. [6l an increase in the particle charge e < 
corresponds to an increase in the velocity v'^ . As the orbital radius decreases, the velocity 
increases until it reaches the limiting value = 1 which corresponds to the velocity of a 
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FIG. 5: The positive solution of the hnear velocity is plotted as a function of the radial distance r/M 
for the parameter choice Q/M = 0.6 and e = — 3 so that /M w 2.74 and the outer horizon is located at 
r+/M = 1.8. The geodesic velocity Vg is also shown (gray curve). The shaded region (r < r+) is forbidden. 
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FIG. 6: The positive solution of the linear velocity is plotted as a function of the radial distance r/M for 
the parameter choice Q/M ~ 0.6 and different values of e = —5 (black curve), e ~ —3 (thick black curve), 
and e = —0.5 (dashed curve). The geodesic velocity Vg for e = is also shown (gray curve). The choice of 
parameters implies that r^ /M k, 2.74 and the outer horizon is located at r^/M = 1.8. The shaded region 
is forbidden. For r > r+ it holds that vf > Vg. 
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FIG. 7: The energy E/ ^ and angular momentum L* = L/{ijlM) of a charged particle of charge-to-mass 
ratio e moving in the field of a RN black hole with charge Q and mass M are plotted as functions of the 
radial distance r/M for the parameter choice Q/M = 0.6 and e = —3, with /M k, 2.7 A and the outer 
horizon located at r+/M = 1.8. The shaded region is forbidden. 

photon. This fact can be seen also in Fig. [TJ where the energy and angular momentum 
for circular orbits are plotted in terms of the distance r. Clearly, this graphic shows that to 
reach the photon orbit at r = r+, the particles must acquire and infinity amount of energy 
and angular momentum. In Fig. [8] we analyze the behavior of the particle's energy and 
angular momentum in terms of the specific charge e. It follows that both quantities decrease 
as the value of lei decreases. 




2. Case e = 



The solutions are the geodesic velocities z/ = ±z/g in the range r > r^. This case has 



been studied in detail in (l2 |. 



3. Case e > 



Depending on the explicit values of the parameters Q and e and the radial coordinate r, 
it is necessary to analyze several subcases. 

a) e < M/Q and r > r^. 
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FIG. 8: The energy E/^ and angular momentum L* = L/{fiM) of a charged particle of charge-to-mass 
ratio e moving in the field of a RN black hole with charge Q and mass M are plotted as functions of the 
radial distance r/M for the parameter choice Q/AI0.6 and e = —3 (solid curves), e = —5 (dashed curves), 
e — —0.5 (dotted curves). Here /M 2.74 and the outer horizon is located at r+/A/ = 1.8. The shaded 
region is forbidden. The energy and angular momentum decrease as |e| decreases. 

There are two different branches for both signs of the hnear velocity: u = ±z/^ in the 
range ri < r < , and u = in the whole range r > ri. The two branches join at 
r = ri, where = v~ = z/^a/A, as shown in Fig. [9l First we note that in this case for 
r > it always holds that i'~ < Ug. This means that, at fixed orbital radius, charged 
test particles possess a smaller orbital velocity than that of neutral test particles in 
the same orbit. This is in accordance to the fact that in this case, a black hole with 
eQ > 0, the attractive gravitational force is balanced by the repulsive electromagnetic 
force. In the region r > r^, the orbital velocity increases as the radius approaches 
the value (see Fig. [9]). The interval ri < r < presents a much more complex 
dynamical structure. First we note that, due to the Coulomb repulsive force, charged 
particle orbits are allowed in a region which is forbidden for neutral test particles. 
This is an interesting result leading to the possibility of accretion disks in which 
the innermost part forms a ring of charged particles only. Indeed, suppose that an 
accretion disk around a RN black hole is made of neutral and charged test particles. 
Then, the accretion disk can exist only in the region r > r/ with a ring of charged 
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FIG. 9: The positive solution of the linear velocity uf- is plotted as a function of the radial distance r/M 
in the region [1.8,5] (left graphic) and [2.5,3] (right graphic). Here Q/AI = 0.6 and e = 1.2, so that 
ri/M = 2.68 and r+ = 2.737A/. For the chosen parameters we have that e = 3.25 and e/ = 4.12. The region 
within the interval [n,r+] is forbidden for neutral particles. 




2.6 2.8 3.0 3.2 3.4 3.6 2.68 2.69 2.70 2.71 2.72 2.73 2.74 

r/M r/M 



FIG. 10: The energy E/ ^ and angular momentum L* = L/{^M) of a charged particle of charge-to-mass 
ratio e moving along circular orbits in a Reissner-Nordstrom black hole of charge Q and mass M are plotted 
in terms of the radial distance r/M in the range [2.6,3.8] (left graphic) [2.68,2.74] (right graphic). Here 
Q/M ~ 0.6 and e = 1.2, so that ri/M = 2.68 and r^ /M = 2.737. For the chosen parameters we have that 
£ = 3.25 e; = 4.12. The shaded region is forbidden for any particles. 

particles in the interval [r/,r+). Outside the exterior radius of the ring (r > r+), 
the disk can be composed of neutral and charged particles. This situation can also 
be read from Fig. [10] where the energy and the angular momentum are plotted as 
functions of the radial distance r/M. 



b) M /Q < e <e and ri < r < r^. 
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FIG. 11: Radius ~ ?'+ (black curve) and ri — (gray curve), are plotted as function of e for Q = 0.6Af. 
r+ = for e = e w 3.25. 
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FIG. 12: Left graphic: The positive solution of the linear velocity v is plotted as a function of the radial 
distance r/M. Right graphic: The energy E/ ^ and angular momentum L* = L/{iiM) of a charged particle 
of charge-to-mass ratio e are plotted in terms of r/M. The parameter choice is Q/AI = 0.6 and e = 2.1. 
Then, ri/M = 2.56, r+ /M = 2.737, and rjM = 3.99. Moreover, for this choice e = 3.25 and ti = 4.12. The 
shaded region is forbidden. 

Since r < r^, one has that e/eo < 1, implying that both solutions v'^ and v~ can exist. 
There are two different branches for both signs: v = in the range ri < r < r+, 
and in the entire range ri < r < Ts- The two branches join at r = r;. Note that 
for increasing values of e, the radius decreases and approaches r^, reaching it at 
e = e, and as e tends to infinity tends to the outer horizon r+ (see Fig. [TT]) . In 
particular, the interaction between the attractive gravitational force and the Coulomb 
force generates a zone r/ < r < r+ in which only charged test particles can move 
along circular trajectories while neutral particles are allowed in the region r > r (see 
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FIG. 13: Left graphic: The positive solution of the hnear velocity i' is plotted as a function of the radial 
distance r/M. Right graphic: The energy E/fi and angular momentum L* = L/{iiM) of a charged particle 
of charge-to-mass ratio e ~ 3.8 moving in a RN spacetime with Q/AI ~ 0.6 are plotted in terms of the radial 
distance r/M. For this choice of parameters the radii are ri/M — 1.98, r^/M — 2.11, and r^/M — 2.737 
whereas the charge parameters are e = 3.25 and ej = 4.12. 

Fig. [12]). This result again could be used to construct around black holes accretion 
disks with rings made of charged particles. 



c) e < e < ei and < r < r+. 

The solution z/~ for the linear velocity is not allowed whereas the solution z/+ is valid 
in the entire range. In fact, the condition r > implies that e/eo > 1, and therefore 
— 1 + (e/eo)^ > A^, so that the condition (138|) for the existence of velocities is 
satisfied for the plus sign only. Therefore, the solutions are given by = ±1/^ in the 
entire range as shown in Fig. [131 At the radius orbit r = r^, the angular momentum 
and the velocity of the test particle vanish, indicating that the particle remains at 
rest with respect to static observers located at infinity. In the region < r < r+ only 
charged particles can move along circular trajectories. 

d) e > ei and < r < r+. 

In this case the radius ri does not exist. The solutions are the velocities u = ±z^^ in 
the entire range. Note that for e — !■ oo one has that — > r+. Also in this case we note 
that neutral particles can stay in circular orbits with a velocity Ug only in the region 
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FIG. 14: Left graphic: The positive solution of the hnear velocity i/ is plotted as a function of the radial 
distance r/M. Right graphic: The energy E/fi and angular momentum L* = L/{^M) of a charged particle 
of charge-to-mass ratio e = 7, moving in the field of a RN black hole with Q/M = 0.6, are plotted in terms 
of the radial distance r/M. For this parameter choice Vg/M — 1.88, e = 3.25, and — 4.12. 

r > r+ whereas charged test particles are allowed within the interval < r < r+, as 
shown in Fig. [TH Clearly, for charged and neutral test particles the circular orbit at 
r = r+ corresponds to a limiting orbit. 

A. Stability 

To analyze the stability of circular orbits for charged test particles in a RN black hole we 
must consider the condition f fT6|) which leads to the Eqs. f|T7|) . f|T71) . and f|T8|) . So the stability 
of circular orbits strongly depends on the sign of {eQ). The case eQ < is illustrated in Fig. 
[TSl where the radius of the last stable circular orbit ri^co is plotted for two different values 
of e as a function of Q/M . It can be seen that the energy and angular momentum of the 
particles decrease as the value of Q/M increases. These graphics also include the radius of 
the outer horizon r+ and the radius r+ which determines the last (unstable) circular orbit of 
neutral particles. In Sec. Illlt we found that circular orbits for charged particles are allowed 
also inside the radius r+ for certain values of the parameters; however, since r+ < ris^o, we 
conclude that all those orbits must be unstable. From Fig. US] we see that for Q = and 
e = 0, the well-known result for the Schwarzschild case, ri^co = QM, is recovered. Also in 
the limiting case Q = M and e = 0, we recover the value of Vi^^^ = 4M for neutral particles 
moving along circular orbits in an extreme BN black hole. In general, as the value of |e| 
increases we see that the value of ri^^^ increases as well. This behavior resembles the case 
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FIG. 15: The radius of the last stable circular orbit Hsco in a RN black hole of mass M and charge Q for a 
particle with ratio e ~ —0.2 (left plot) and e ~ —1.5 (right plot). Numbers close to the point represent the 
energy E/ ^ of the last stable circular orbits at that point. Underlined numbers represent the corresponding 
angular momentum L/{Mfi). Stable orbits are possible only for r > Visco- For comparison we also include 
the curves for the radii r+ and rt . 
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FIG. 16: The radius risco of the last stable circular orbit in a RN black with charge-to-mass ratio Q/M for 
selected values of the charge-to-mass ratio e of the test particle. Only the case eQ < is illustrated. Stable 
orbits are possible only for r > ri^co 



of the radius of the last stable orbit for neutral test particles |12l . Il3l | . Indeed, in the case 
eQ < the attractive Coulomb force reinforces the attractive gravitational force so that the 
general structure remains unchanged. We also can expect that an increase in the charge of 
the particle |e| produces an increase in the velocity of the stable circular orbits. In fact, this 
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FIG. 17: The radius of the last stable circular orbit risco (solid curve) for a charged test particle with 
e = 7, in a RN black hole with charge Q and mass Af , is plotted as a function of the ratio Q/M . Other 



curves arc the outer horizon radius ?■+ = M + \J — and the radii ?-+ = [3A/ + ^ (9Af^ — 8Q^)]/2, 
= [ca/M^— Q2^^2— T + il/(e2 - 1)1, = ^ + i v^SA^^ - gg^ - Q2g2_ Shaded and dark 

regions are forbidden for timclike particles. Stable orbits are possible only for r > risco- 

can be seen explicitly from Eq. (130|) and Fig. |6l It then follows that the energy and angular 
momentum of the charged test particle increases as the value of |e| increases. 

The case of eQ > is illustrated in Figs. [17] and [181 The situation is very different 
from the case of neutral particles or charged particles with eQ < 0. Indeed, in this case 
the Coulomb force is repulsive and leads to a non trivial interaction with the attractive 



gravitational force, see also [32|-[4l|. It is necessary to analyze two different subcases. The 
first subcase for e > 1 is illustrated in Fig. [T71 while the second one for < e < 1 is depicted 
in Fig. [TH We can see that in the case < e < 1 the stability regions are similar to 
those found in the case e < (cf. Figs. [TSl and [TSl) . This means that for weakly-charged 
test particles, < e < 1, it always exists a stable circular orbit and Vi^^^ > 4M, where the 
equality holds for an extreme black hole. On the contrary, in the case e > 1 there are regions 
of Q and e in which stable circular orbits cannot exist at all. As can be seen from Fig. [T7[ 



charged particles moving along circular orbits with radii located within the region r < r+ 
or r < r<. must be unstable. 

We conclude that the ring structure of the hypothetical accretion disks around a RN 
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FIG. 18: The radius of the last stable chxular orbit ri^co (solid curve) for a charged test particle with 
e = 0.5, in a RN black hole with charge Q and mass M, is plotted as a function of the ratio Q/M. Other 
curves are the outer horizon radius and the radius . Shaded and dark regions are forbidden for timelike 
particles. Stable orbits are possible only for r > ri^co- 

black hole mentioned in Sec. IIIII must be unstable. 



IV. NAKED SINGULARITIES 

The effective potential V± given in Eq. (9) in the case of naked singularities (M^ < Q"^) 
is plotted in Figs. ( IT9] - [22|) in terms of the radial coordinate r/M for selected values of the 



ratio e and the angular momentum L/{Mfi) of the test particle, see also 25|, |3l|, |42N47|. 
The effective potential profile strongly depends on the sign of eQ. Moreover, the cases with 
|e| < 1 and with |e| > 1 must be explored separately. 

Fig. [20] shows the effective potential for a particle of charge-to-mass e in the range 
[—10,-1]. The presence of minima (stable circular orbits) in the effective potential with 
negative energy states is evident. Moreover, we note that the minimum of each potential 
decreases as |e| increases. This fact is due to the attractive and repulsive effects of the 



gravitational and electric forces (32l-l4l| . In Fig. |2T] the effective potential is plotted for 
negative and positive values of the charge-to-mass ratio e. We see that for a fixed value of 
the radial coordinate and the angular momentum of the particle, the value of the potential 
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FIG. 19: The effective potential for a charged particle with charge-to-mass ratio e in a RN naked singularity 
of charge Q and mass M is plotted as a function of the radius r jM for fixed values of the angular momentum 
L* = L/{iiM). Black curves represent the positive solution V+ while gray curves correspond to V-. The 
boldfaced points denote the minima of the potentials. In upper left plot, the parameter choice is Q/M = 2 
and e = —1.5; the upper right plot is for Q/M = 2 and e = — 5 while the bottom plot corresponds to the 
choice Q/M = 1.5 and e = -0.2. 

V increases as the value of e increases. In the Fig. [22] we plot the effective potential for a 
fixed Q/M as function of the radial coordinate and the angular momentum for two different 
cases, e = 0.1 and e = —0.1. We can see that in the first case the presence of a repulsive 
Coulomb force reduces the value of the radius of the last stable circular orbit for a fixed 
angular momentum. We note the existence of stable "circular" orbits with L = at which 
the particle is at rest with respect to static observers located at infinity. 

Negative energy states are possible only in the case eQ < 0. The region in which the 
solution V+ has negative energy states is 

< r < M + ^/A'P - (1 _ e2) fo^ e < -1 , (49) 
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FIG. 20: The effective potential of a RN naked singularity with Q/M = 2 for a particle with charge-to-mass 
ratio e in the range [—10, —1] and angular momentum L/(i\//i) = 4. 
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FIG. 21: The effective potential of a RN naked singularity with Q/M = 3/2 for a particle with charge- 
to-mass ratio e in the range [—2, +2] and angular momentum L/{Mfi) = 4. 



and 



< r < rl' for < L < Lg , e < -1 



Ti <r <rl for < L < Lq , -1 < e < -W 1 



M2 



where 



- 1 . 



(50) 
(51) 

(52) 



In general, for a particle in circular motion with radius tq and charge-to-mass ratio 
e, around a RN naked singularity with charge Q and mass M, the corresponding angular 
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FIG. 22: The effective potential V+ of a RN naked singularity with Q/M = 3/2 for a charged particle is 
plotted for different values of the angular momentum L* = L/{M^). The left plot corresponds to the ratio 
e = 0.1 while the right one is for e = —0.1. For e = 0.1 there is a minimum, Vmin ~ 0.81, at r,„i„ « 2.52A/ for 
L* = 0, and a minimum, Vrrim « 0.96, at « 29Af for L* = 5. For e = —0.1 the minimum, Vmin ~ 0.67, 



is located at 



2.02A'/ for L* = 0, and at 



20.8Af with Vrr, 



0.97 for L* = 5. 



momentum must be chosen as 



,2 <''0 



(53) 



11^ - \rl - 2Mro + 

in order for negative energy states to exist. 

The conditions for circular motion around a RN naked singularity are determined by 
Eq. ( |TOl) which can be used to find the energy and angular momentum of the test particle. 
Indeed, Eqs. fllll) and (IT^ define the angular momentum and the energy E^, respectively, 
in terms of r/M, Q/M, and e. The explicit dependence of these parameters makes it 
necessary to investigate several intervals of values. To this end, it is useful to introduce the 
following notation 



3M 



r, = 



2 ±^V9M2-8g2_Q2e2 



y/2Q 



and 



y/2Q 



3Af2 - 2g2 ± Af ^9Af2 _ 8g2 



(54) 
(55) 



We note that 



limr^ 

e->-0 ^ 



M 



(56) 



(57) 
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FIG. 23: The charge parameters e; (black sohd curve), e_ (gray sohd curve), e+ (dashed curve), e_ (dotted 
curve), and e+ (dotdashed curve) as functions of the charge-to-mass ratio of the RN naked singularity. The 
special lines Q/M = 5/ (2^6) w 1.02, Q/M = 3^6/7 w 1.05, and Q/M = i/o/S w 1.06 are also plotted. 



which corresponds to the classical radius of a mass M with charge Q, see for example 48l. |49| . 
and 



lim r," 



3M 



2 ±iv9M2-8g2 



(58) 



which represents the limiting radius at which neutral particles can be in circular motion 
around a RN naked singularity 12|. 

The behavior of the charge parameters defined above is depicted in Fig. |23] in terms of 
the ratio Q/M > 1. It follows from Fig. |23]that it is necessary to consider the following 
intervals: 



Q/M G (1,5/(276)], 

Q/M G (5/(2v^),(3V6)/7], 

Q/M G i{3V6)/7, ^9/8], 

Q/M G [v^,oo). 



(59) 
(60) 
(61) 
(62) 



Our approach consists in analyzing the conditions for the existence of circular orbits by 
using the expressions for the angular momentum, Eq. f lTTj) . and the energy, Eq. f|T2|) . of the 
particle together with the expressions for the velocity obtained in Sec. Illli We consider 
separately the case e > in Sees. IIV Al and lIVBt and e < in Sees. IIV CI and IIVDI In 
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FIG. 24: The case e > 1. The energy (black curve) and angular momentum (gray curve) for a test particle 
with charge-to-mass ratio e = 2 in a RN naked singularity with Q ~ 1.06Af . Circular orbits exist in the 
interval < r < r+, where r~ = 1.04196Af and r+ = 1.95804il/. 

the Appendix |X] we present equivalent results by using the alternative method of the proper 
linear velocity of test particles in an orthonormal frame as formulated in Sec. Illli 

A. Case e > 1 

For e > the condition (l33ll implies in general that r > = jM. Imposing this 
constraint on Eqs.( |TT]) and (fT2l) . we obtain the following results for timelike orbits. For 
e > 1 and M < Q < a/9/8M circular orbits exist with angular momentum L = L"*" in the 
interval r~ < r < r+, while for Q > a/9/8 Af no circular orbits exist (see Fig. [2^ . Clearly, 
the energy and angular momentum of circular orbits diverge as r approaches the limiting 
orbits at r^. This means that charged test particles located in the region r~ < r < need 
to acquire an infinite amount of energy to reach the orbits at r^. The energy of the states 
is always positive. A hypothetical accretion disk would consist in this case of a charged 
ring of inner radius r~ and outer radius r+, surrounded by a disk of neutral particles. The 
boundary r = r;^ in this case would be a lightlike hypersurface. 

Since for eQ > the Coulomb interaction is repulsive, the situation characterized by 
the values for Q > ^9/8M and e > 1 corresponds to a repulsive electromagnetic effect 
that cannot be balanced by the attractive gravitational interaction. We note that the case 
Q > a/9/8M and e > 1 could be associated to the realistic configuration of a positive ion 
or a positron in the background of a RN naked singularity. 
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B. Case < e < 1 



It turns out that in this case it is necessary to consider separately each of the four different 
regions for the ratio Q/M that follow from Fig. [221 Moreover, in each region of Q/M it 
is also necessary to consider the value of e for each of the zones determined by the charge 
parameters e/, e±, and e±, as shown in Fig. [231 We analyzed all the resulting cases in detail 
and found the values of the energy and angular momentum of charged test particles in all 
the intervals where circular motion is allowed. We summarize the results as follows. 

There is always a minimum radius rmin at which circular motion is allowed. We found 
that either Tmin = oi' ^mm = ■ Usually, at the radius r+ the test particle acquires a zero 
angular momentum so that a static observer at infinity would consider the particle as being 
at rest. Furthermore, at the radius r~ the energy of the test particle diverges, indicating 
that the hypersurface r = r~ is lightlike. In the simplest case, circular orbits are allowed 
in the infinite interval [rmin-, oo) so that, at any given radius greater than rmin-, it is always 
possible to have a charged test particle moving on a circular trajectory. Sometimes, inside 
the infinite interval [r^jn, oo), there exists a lightlike hypersurface situated at r+ > r^m- 

Another possible structure is that of a finite region filled with charged particles within the 
spatial interval (r^m = , Vmax = r^)- This region is usually surrounded by an empty finite 
region in which no motion is allowed. Outside the empty region, we find a zone of allowed 
circular motion in which either only neutral particles or neutral and charged particles can 
exist in circular motion. Clearly, this spatial configuration formed by two separated regions 
in which circular motion is allowed, could be used to build with test particles an accretion 
disk of disconnected rings. A particular example of this case is illustrated in Fig. [25] 



C. Case e < -1 



The contribution of the electromagnetic interaction in this case is always attractive. 
Hence, the only repulsive force to balance the attractive effects of the gravitational and 
Coulomb interactions can be generated only by the RN naked singularity. This case therefore 
can be compared with the neutral test particle motion as studied in [12I, 13|. Then, it is 



convenient, as in the case of a neutral test particle, to consider the two regions 
and M < g < ^/o/SM separately. 
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FIG. 25: Case: M < Q < 5/(2\/6)A/ and e+ < e < e;. Parameter choice: Q ^ I.OIM and e = 0.902. 



Then ei = 0.907, ?+ = 0.8963, r+ = 1.44942M, = 1.04196M, r 



r-+ - 



1.95804Af, rr 



1.45192M, and 



r'^ = 1.548077Af . Ch-cular orbits exist with angular momentum L = (gray curves) and energy E = 
(black curves) in r;^; < r < rf (upper left plot); L = in rf < r < r'[ (upper right plot) and r'^ < r < r+ 
(bottom left plot); L = in ?- > r+ (bottom right plot). 



For e < — 1 and for Q > ^9/8M circular orbits with L = always exist for r > 
(in fact, however, one has to consider also the limit r > for the existence of timelike 
trajectories). This case is illustrated in Fig. [261 where the presence of orbits with negative 
energy states is evident. 



For M < Q < a/9/8M circular orbits exist with L = in < r < and r > r+ (see 
Fig. [27|) . We note that for neutral test particles in the region M <Q< ^/9/8M, (stable) 
circular orbits are possible for r > r* = Q'^/M. At r = r*, the angular momentum of the 
particle vanishes 12|. On the contrary, charged test particles with e < — 1 can move along 
circular orbits also in the region (0,r*]. The value of the energy on circular orbits increases 
as r approaches r = 0. However, the angular momentum, as seen by an observer located at 
infinity, decreases as the radius of the orbit decreases. In the region M <Q< ^/9/8M, two 
limiting orbits appear at r^, as in the neutral particle case 12|. 
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FIG. 26: Case: e < — 1 and Q > ■\/9/8Af . Parameter choice: Q = 2M and e = —2. Circular orbits exist 
with angular momentum L = (gray curve) and energy E = (black curve). 
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FIG. 27: Case: e < -1 and Af < Q < yJoJSM. Parameter choice: Q = I.QIM and e = -2. Then, 
r~ = 1.04196 A/ and r+ = 1.95804A/. Circular orbits exist with angular momentum L = (gray curve) 
and energy E = E^ (black curve) inO < r < and r > r+. 

D. Case -1 < e < 

For this range of the ratio e, it is also convenient to analyze separately the two cases 
Q > ^/9/8M and M < Q < ^/9/8M. In each case it is necessary to analyze the explicit 
value of e with respect to the ratio M/ Q. Several cases arise in which we must find the regions 
where circular motion is allowed and the value of the angular momentum and energy of the 
rotating charged test particles. 

We summarize the results in the following manner. There are two different configurations 
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FIG. 28: Case: M < Q < ^/9/8AI and -M/Q < e < 0. Parameter choice: Q = 1.05M and e = -0.2. 
Then = 1.28787^, r+ = 1.71213A/, and = 1.03487Af . Circular orbits exist with angular momentum 
L = (gray curve) and energy E = (black curve) in < r < (left plot) and in r > r+ (right 
plot). For r = , L = 0. 

for the regions in which circular motion of charged test particles is allowed. The first one 
arises in the case Q > a/9/8M, and consists in a continuous region that extends from a 
minimum radius Tmin to infinity, in principle. The explicit value of the minimum radius 
depends on the value of e and can be either r~, r+, or rmin = Q^/{2M). In general, we find 
that particles standing on the minimum radius are characterized by L = 0, i. e., they are 
static with respect to a non-rotating observer located at infinity. 

The second configuration appears for M < Q < ^9/8M. It also extends from r^m to 
infinity, but inside it there is a forbidden region delimited by the radii r~ and r^. The 
configuration is therefore composed of two disconnected regions. At the minimum radius, 
test particles are characterized by L = 0. On the boundaries (r^) of the interior forbidden 
region only photons can stand on circular orbits. A particular example of this case is 
presented in Fig. [281 

E. Stability 

To explore the stability properties of the circular motion of charged test particles in a RN 
naked singularity, it is necessary to investigate the equation flTBl) or, equivalently, Eqs. fll7l) . 
(fTTll . and (JTSl) . considering the different values for e and Q/M > 1. We can distinguish two 
different cases, |e| > 1 and < |e| < 1. Let us consider the case |e| > 1. In particular, as it 
was shown in Sec. IIV At for e > 1 and M < Q < ^9/8M circular orbits exist with L = L"*" 
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FIG. 29: The radius of the last stable circular orbit ri^co (gray curve) of a charged particle with ratio 
e = +7 in a RN naked singularity with ratio Q/M G [1,1.2]. The radii r, = Q^/M and r = = 
[iM ± ^9A/2 _ 8Q2J/2 are also plotted. Circular orbits exist only in the interval KQ/M < ^/9/8. The 
shaded region is forbidden for timelike particles. Stable orbits are located in the region r > ri^co- 

in the interval r~ < r < whereas no circular orbits exist for e > 1 and Q > a/9/8M. 
For this particular case, a numerical analysis of condition f lT6|l leads to the conclusion that 
a circular orbit is stable only if its radius tq satisfies the condition vq > Vi^^^, where Vi^^^ is 
depicted in Fig. [221 We see that in general the radius of the last stable circular orbit is 
located inside the interval {r~,r^). It then follows that the only stable region is determined 
by the interval ri^^o < r < r+. 

Consider now the case e < — 1. The numerical investigation of the condition f fT6|) for 
the last stable circular orbit shows that in this case there are two solutions r^^^ such that 
f^co — ^ico^ where the equality is valid for Q/M ^ 1.72. Moreover, for Q/M = a/9/8 
we obtain that r~^^ — = ■ This situation is illustrated in Fig. [301 Stable orbits 
corresponds to points located outside the region delimited by the curves r = r+^, r = r~^^, 
and the axis Q/M = 1. On the other hand, we found in Sec. II V CI that for e < —1 and 
1 < Q/M < a/9/8 circular orbits exist in the interval < r < r ~ and r > r^. It then follows 
that the region of stability corresponds in this case to two disconnected zones determined 
by < r < r~ and r > r^^^. Moreover, we established in Sec. IIV CI that for e < — 1 
and ^J9/8 < Q/M circular orbits always exist for r > 0. Consequently, in the interval 
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1.0 1.2 1.4 1.6 1.8 1.00 1.01 1.02 1.03 1.04 1.05 1.06 

Q/M Q/M 

FIG. 30: The radius of the last stable circular orbit r^^^ (black curves) of a charged particle with ratio 
e = —7 in a RN naked singularity with ratio Q/M <E [1,1.8]. The radii r*, and are also plotted for 
comparison. In the shaded region no circular orbits can exist. Stable circular orbits are situated outside the 
region with boundaries r^^^, ''i7coi ^-nd the vertical axis Q/M = 1. 

a/9/8 < Q/M < 1.72, the stable circular orbits are located in the two disconnected regions 
defined by < r < r^^^^ and r > r+^. Finally, for Q/M > 1.72 all the circular orbits are 
stable (see Fig. [30]). 

The case < |e| < 1 is much more complex, and needs to be described for different 
subcases following the classification of orbital regions traced in Sec. IIVBI for the case 
< e < 1, and in Sec. IIVDI for the case — 1 < e < 0. The results for the specific ratio 
e = 0.5 are given in Fig. [31] and for e = —0.5 in Fig. [32] In general, we find that the results 
are similar to those obtained for the case e < — 1. Indeed, the zone of stability consists of 
either one connected region or two disconnected regions. The explicit value of the radii that 
determine the boundaries of the stability regions depend on the particular values of the ratio 
Q/M. 
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FIG. 31: The radius of the last stable circular orbit r^^^ (black and grey curves) of a charged particle 



with ratio e = 0.5 in a RN naked singularity with ratio Q/M <S [1,1.1]. Also plotted: 



[3M ± 



+ — 



v/e2(e2_i)(M2-Q2)^/(g2 _ 



, r± = ± i^9M2-8Q2_Q2g2^ 



and r, = Q^/M. Regions of stability are: for Q > ^/9/8M in r > r^, for (3^6/7)1/ < Q < ^/9/8M 
exist stable orbits in r~ < r, for (5/(2-\/6))A/ < Q < (3\/6/7)A/ exist stable orbits in < r. For 



M <Q < (5/(2^/6))M stable orbits are located in r > 



Isco ' 
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1.00 
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1.10 
Q/M 



1.15 



1.20 



FIG. 32: The radius of the last stable circular orbit r^^^ (black and grey curves) of a charged particle 
with ratio e = —0.5 in a RN naked sing ularity with ratio Q/M e [1,1.2]. Also plotted: r± = [3A/ ± 



^)]/2, r^ 



3M 
2 



± iV9^ 



2,2 



^e2(e2_l)(M2-Q2)^./(g2 _ ^-^ 

and r, = Q^/M. Shaded regions are forbidden. Regions of stability are: for Q > {y^9/8)M stable circular 
orbits exist in rf < r < rf^^^, and r > r^^^^- For M < Q < (■\/9/8)Af stable circular orbits exist in 



rf < r < , and r > r+. 
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V. CONCLUSIONS 



In this work, we explored the motion of charged test particles along circular orbits in the 
spacetime described by the Reissner-Nordstrom (RN) metric. We performed a very detailed 
discussion of all the regions of the spacetime where circular orbits are allowed, using as 
parameters the charge-to-mass ratio Q/M of the source of gravity and the charge-to-mass 
ratio e = g/yU of the test particle. Depending on the value of Q/M, two major cases must 
be considered: The black hole case, \Q/M\ < 1, and the naked singularity case, \Q/M\ > 1. 
Moreover, we found out that the two cases |e| < 1 and |e| > 1 must also be investigated 
separately. Whereas the investigation of the motion of charged test particles with |e| > 1 can 
be carried out in a relatively simple manner, the case with |e| < 1 is much more complex, 
because it is necessary to consider various subcases which depend on the explicit value of e 
in this interval. 

To perform the analysis of circular motion of charged test particles in this gravitational 
field we use two different methods. The first one consists in using constants of motion to 
reduce the equations of motion to a single first-order differential equation for a particle 
moving in an effective potential. The properties of this effective potential are then used to 
find the conditions under which circular motion is possible. The second approach uses a 
local orthonormal frame to introduce a "local proper linear velocity" for the test particle. 
The conditions for this velocity to be timelike are then used to determine the regions of 
space where circular orbits are allowed. The results of both methods are equivalent and, 
in fact, for the sake of simplicity it is sometimes convenient to use a combination of both 
approaches. In this work, we analyzed in detail the conditions for the existence of circular 
orbits and found all the solutions for all the regions of space in the case of black holes and 
naked singularities. 

To formulate the main results of this work in a plausible manner, let us suppose that 
an accretion disk around a RN gravitational source can be made of test particles moving 
along circular orbits 50||. Then, in the case of black hole we find two different types of 
accretion disks made of charged test particles. The first type consists of a disk that begins 
at a minimum radius R and can extend to infinity, in principle. In the second possible 
configuration, we find a circular ring of charged particles with radii {rint,rext), surrounded 
by the disk, i. e., with r^rct < R- For certain choices of the parameter e the exterior disk 
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might be composed only of neutral particles. A study of the stability of circular orbits 
shows that the second structure of a ring plus a disk is highly unstable. This means that 
test particles in stable circular motion around RN black holes can be put together to form 
only a single disk that can, in principle, extend to infinity. 

In the case of RN naked singularities we find the same two types of accretion disks. 
The explicit values of the radii Tmin, ^ext, and R depend on the values of the ratios e and 
Q/M, and differ significantly from the case of black holes. In fact, we find that the case of 
naked singularities offers a much richer combination of values of the charge-to-mass ratios 
for which it is possible to find a structure composed of an interior ring plus an exterior 
disk. A study of the stability of this specific situation shows that for certain quite general 
combinations of the parameters the configuration is stable. This result implies that test 
particles in stable circular motion around RN naked singularities can be put together to 
form either a single disk that can extend, in principle, to infinity or a configuration of an 
interior ring with an exterior disk. This is the main difference between black holes and naked 
singularities from the viewpoint of these hypothetical accretion disks made of test particles. 

The question arises whether it is possible to generalize these results to the case of more 
realistic accretion disks around more general gravitational sources, taking into account, for 



instance, the rotation of the central body, 5l|, |52|. It seems reasonable to expect that in 



the case of Kerr and Kerr-Newman naked singularities, regions can be found where stable 
circular motion is not allowed so that an accretion disk around such an object would exhibit 
a discontinuous structure. Indeed, some preliminary calculations of circular geodesies in 
the field of rotating compact objects support this expectation. Thus, we can conjecture 
that the discontinuities in the accretion disks around naked singularities are a consequence 
of the intensity of the repulsive gravity effects that characterize these speculative objects. 
Furthermore, it was recently proposed that static compact objects with quadrupole moment 
can be interpreted as describing the exterior gravitational field of naked singularities 531.154 1. 
It would be interesting to test the above conjecture in this relatively simple case in which 
rotation is absent. If the conjecture turns out to be true, it would give us the possibility 
of distinguishing between black holes and naked singularities by observing their accretion 
disks. 
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Appendix A: Velocity of test particles in a RN naked singularity 

In this Appendix we explore charged test particles in circular motion in a RN naked sin- 
gularity by using the tetrad formalism, as developed in Sec. IIIII for the black hole analysis. 
In Sec. llVt we studied the timelike circular motion in the naked singularity case by analyz- 
ing directly the existence conditions for the energy, Eq. f fT2|) . and the angular momentum, 
Eq. (ITT]) . Here we use the formalism of "local proper linear velocity" as measured by an 
observer attached to an orthonormal frame. The results are equivalent to those obtained by 
using the expressions for the energy and angular momentum. 

In Sec. mil we showed that the linear velocity of a test particle in a RN spacetime can 
be written as 

± 



1/2 



A ± - 1 + (e/eo)2j , (Al) 
where 

Then, the conditions for the existence of timelike velocities are 



A^ - 1 + (e/eo)2 > , (A3) 
A± VA2-l + (e/eo)2 > , (A4) 
{utr < 1 . (A5) 

We first note that, in the case of a naked singularity, these conditions can be satisfied only 
for r > QVM. 

For e > 1 and e < — 1 the solutions are the geodesic velocities u = it'^. In fact, in this 
case, condition (1A4I) with the minus sign is no more satisfied. On the other hand, conditions 
f lA3|) . (]A4p . and (lASP imply that circular timelike orbits exist ior Q/M > \/9/8 in the entire 



range r > Q'^/M. For 1 < Q/M < \/9/8 circular orbits are possible in r > Q'^/M and 
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r ^ r± = [3M ± ^9M2^^8Q2]/2. Finally, for Q/M = ^/9/8 timelike circular orbits exist 
for all r > Q'^/M, except at r = (3/2)M. Moreover, the radii r = correspond to photon 
orbits in the RN spacetime (see Fig. |33|) . 

Consider now the case |e| < 1. It is useful to introduce here the following notations: 



e± = -^y^5M2 ± 4Q2 + ^/25M2 - 24g2, (AT) 



and 



TT ^ -iQ^zrW^ [M{e' - 1) ± V62(e2-l)(M2-g2)] . (A8) 

First, consider the case < e < 1. For e > condition ( 133|) implies that r > jM. 
Applying this constraint on conditions flA3l) and flA4l) . we obtain the following results for 
timelike geodesies. 

1. For 1 < Q/Af < 5/(2v/6) the following sub cases occur: 



a) < e < e_ : Fig. [3li 

The velocity v = ±z/^ exists in the range Q'^/M < r < rf and r > r^' with 
r 7^ r^, V = ±z/~ exists in the range < r < r^" and r > r^"^. 



b) e_ < e < e+ : Fig. 

The velocity = ±z^^ exists in the range Q'^/M < r < rf and r > r^' with 
r 7^ r^, z/ = ±z/~ exists in the range r > r^. 

c) e+ < e < Q, : Fig. El 

The velocity u = ±z^^ exists in the range Q'^/M < r < and r > r^' with 
r 7^ r^, z/ = ±z/~ exists in the range r+ < r < r^" and r > r'^. 

d) Q < e < 1 : Fig. EH 

The solutions are the geodesic velocities u = in the range r > Q'^/M with 
r 7^ r^. The solution u = ±z/~ exists for ei < e < M/Q in the range r > r+. 
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2. For 5/(2^6) < Q/M < the following subcases occur: 

a) < e < : Fig. Elo 

The velocity v = ±1/^ exists in the range Q'^/M < r < and r > r'l' with 
r 7^ r^, V = ±z/~ exists in the range < r < rf and r > r^. 

b) ez < e < 1 : Fig. IMJi 

The velocity z/ = exists in the range r > Q'^/M, v = ±z/~ exists in the range 
r > r+. 

3. Q/M > ^/9/8 : Figs. EE] and [37] 

The velocity u = exists in the range r > Q'^/M for Q/M > a/9/8 whereas for 
Q/M = y/g/S this is a solution in r/M > 9/8 with r/M ^ 3/2, u = exists for 
< e < M/Q in the range r > rf. 

The results for — 1 < e < are summarized below. 

1. For 1 < Q/M < 5/(276) the following sub cases occur: 

a) For — 1 < e < — e;, the velocity u = ±uf exists in the range r > Q'^/M with 

r 7^ r^, = ±i/~ exists for —{M/Q) < e < — e; in the range r > (see Fig. 

Ell). 

b) For —ei < e < — e+, the solution is z/ = ±i^^ in the range /M < r < rf and 

r > r'f with r 7^ r^, z/ = iz/^ exists in the range rf < r < rf and r > rf (see 

Fig. EHb). 

c) For — e+ < e < — e_, the velocity u = ±z/~ exists in the range r > rf. v = ±z/+ 

exists for — e+ < e < — e_ in the range {Q'^/M) < r < rf, and r > r/' with 
r ^ rf, and for e = — the velocity z/+ exists for Q^/M < r < rf and 
r > r^"^ with r 7^ r^. Finally, for Q = 5/(2a/6)M and e = — e+, z/+ exists for 
(Q^/M) <r <rf, and r> rf (see Fig. EHt). 
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FIG. 33: The positive solution of the linear velocity f + is plotted as a function of the radial distance r/M 
for different values of the ratios Q/M and e. The geodesic velocity lyg is also shown (dashed curve). Shaded 
region is forbidden. In (a) the parameter choice is Q/M = 3 and e = 2, with = Q'^/M = 9M. In 
(6) the parameter choice is Q/M = 2 and e = 3, with r, = Q^/M = AM. In (c) the parameter choice 
is Q/M = 1.04 and e = 2, with r* = Q^/M « 1.08A/, r+ = [UI + ^/9M^ - 8Q^]/2 w 1.79A/, and 



[3A/ - y/QAP - 8Q2J/2 w 1.201A/. In (d) the parameter choice is Q/M = ^/9/8 and e = 2, with 
Q^M ^ {9/8)M, r± = [3M ± ^9^/2 - 8Q^]/2 = {3/2)M. 



d) For — e_ < e < 0, the solutions are the geodesic velocities u = ±1^^ in the range 
(Q'^/M) < r < r/^ and r > r^*^ with r ^ r^. The solution = ±u~ exists in 
r+ < r < rf and r > r^"^ (see Fig. 



2. For 5/(2a/6) < Q/M < a/9/8 the following subcases occur: 

a) For — 1 < e < — e;, the velocity z/ = ±1/^ exists in the range r > Q'^/M with 
r 7^ r^, = ±i^~ exists for —{M/Q) < e < — in the range r > (see Fig. 
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(c) (d) 
FIG. 34: The positive solution of the linear velocity is plotted as a function of the radial distance 
r/M for Q/M = 1.01 and different values of the ratio e. In this case rj^ ~ 1.96M, — 1.042iV/ with 
= /AI — 1.02il/, e_ = 0.31, e+ = 0.9, e; w 0.91. The geodesic velocity i^g is also shown (dashed curve). 
Shaded region is forbidden. In (a) the parameter choice is e = 0.2 with rf = 1.05Af, and r'^ = 1.95Af, 
r," = 1.05M. In (6) the parameter choice is e = 0.5. Here r+ = l.llAf, and r,+ = 1.88M, = 1.12M. In 
(c) the parameter choice is e = 0.9. Here rf — l.llM, and r^' — 1.88M, = 1.12M. In {d) the parameter 
choice is e = 0.95. Here r+ = 1.79A/, and rf^ do not exist. 



b) For —ei < e < 0, the velocity u = exists in the range Q'^/M < r < Vi and 
r > r^, r ^ r^, v = ±u~ exists in the range r+ < r < and r > (see Fig. 



3. For Q/M > \/9/8 the velocity u = exists for —{M/Q) < e < in the range 
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(a) 



(b) 



FIG. 35: The positive solution of the linear velocity is plotted as a function of the radial distance 
r/M for Q/M = 1.05 and different values of the ratio e. In this case rj^ = 1.71A/. r.:^ ~ 1.29il/ with 
= Q^/M w 1.102il/, ei w 0.40. The geodesic velocity I'g is also shown (dashed curve). Shaded region is 
forbidden. In (a) the parameter choice is e = 0.7. Here rf = 1.61il/, and rf arc not defined. In (6) the 



parameter choice is e = 0.2. Here rj = 1.2M, and r, = 1.68M, 



1.32M. 





12 3 4 

r/M 

(a) (b) 
FIG. 36: The positive solution of the linear velocity j/^ is plotted as a function of the radial distance r/M 
for Q/M ~ \/9/8 and different values of the ratio e. In this case r+ = r.:^ = 3/2M with r* = 9/8A/, e; = 0. 
The geodesic velocity Vg is also shown (dashed curve). Shaded region is forbidden. In (a) the parameter 
choice is e = 0.2 with rf ~ 1.2M. In (&) the parameter choice is e = 0.7 with rf = 1.72A'/. 



r > r~ 



V = ±i^+ is a solution for Q/M > y^9j8 and -1 < e < in r > Q'^/M 
whereas for Q/M = a/9/8 this is a solution in r/M > 9/8 with r/M ^ 3/2 (see 
Figs. HQ] and E]). 
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FIG. 37: The positive solution of the hnear velocity is plotted as a function of the radial distance r/M 
for Q/M = 2 and different values of the ratio e. In this case r* = Q^/M w AM. The geodesic velocity 
is also shown (dashed curve). Shaded region is forbidden. In (a) the parameter choice is e = 0.7 with 
rf — 1.48M. In (6) the parameter choice is e = 0.2 with rf = 6.2M. 
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(c) 



(d) 



FIG. 38: The positive solution of the linear velocity i/^ is plotted as a function of the radial distance r/M for 
Q/M = 1.01 and different values of the ratio e. In this case r+ = 1.96M, = 1.042M = Q'^/M ^ 1.02M, 
e_ = 0.31, £-1- ~ 0.9, and e; w 0.91. The geodesic velocity i/g is also shown (dashed curve). Shaded region 
is forbidden. In (a) the parameter choice is e = —0.95. Here r+ — 1.79M, and rf^ do not exist. In (6) the 
parameter choice is e = —0.9. Here r+ = l.llM, r'^ = 1.88M, and = 1.12i\/. In (c) the parameter 



choice is e = -0.5. Here r+ = l.llA/, = 1.88Af, and 



1.12A/. In [d) the parameter choice is 



-0.2. Here r+ = 1.05M, 



1.95M, and = 1.05M. 
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(a) 



(b) 



FIG. 39: The positive solution of the hnear velocity is plotted as a function of the radial distance r/M 
for Q/M = 1.05 and different values of the ratio e. In this case r+ = 1.71il/, r;:^ = 1.29M, = Q'^/M ss 
1.102A/, and e; w 0.40. The geodesic velocity Vg is also shown (dashed curve). Shaded region is forbidden. 



In (a) the parameter choice is e = —0.2. Here ?'+ = 1.2A/ 



' 'I 



1.68A/, and = 1.327\/. In (&) the 



parameter choice is e = —0.7. Here r+ = 1.61A/, and r, are not defined. 





12 3 4 

r/M 

(a) (b) 

FIG. 40: The positive solution of the linear velocity is plotted as a function of the radial distance r/M 
for Q/M = ^9/8 and different values of the ratio e. In this case = ^ 3/27\/, r» = 9/8A/, and e; = 0. 
The geodesic velocity Vg is also shown (dashed curve). Shaded region is forbidden. In (a) the parameter 
choice is e = —0.7 with r+ = 1.72A/. In (6) the parameter choice is e = —0.2 with r+ = 1.2M. 
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FIG. 41: The positive solution of the hnear velocity is plotted as a function of the radial distance r/M 
for Q/M = 2 and different values of the ratio e. In this case r» = Q"^ /M ~ 4A/. The geodesic velocity Vg 
is also shown (dashed curve). Shaded region is forbidden. In (a) the parameter choice is e = —0.2. Here 
r+ = 6.2Af. In (6) the parameter choice is e = —0.7. Here r+ = 1.48Af . 
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